Conductance fluctuations in the presence of spin scattering 



m 
O 
o 

(N 

Oh' 
(N 



Jun-ichiro Ohe, Masayuki Yamamoto, and Tomi Ohtsuki 
Department of Physics, Sophia University, Kioi-cho 7-1, Chiyoda-ku, Tokyo 102-8554, Japan 

(Dated: February 2, 2008) 

Electron transport through disordered systems that include spin scatterers is studied numerically. 
We consider three kinds of magnetic impurities: the Ising, the XY and the Heisenberg. By extending 
the transfer matrix method to include the spin degree of freedom, the two terminal conductance 
is calculated. The variance of conductance is halved as the number of spin components of the 
magnetic impurities increases. Application of the Zeeman field in the lead causes a further halving 
of the variance under certain conditions. 

PACS numbers: 73.23.-b, 72.25.Rb, 73.63. Nm 
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I. INTRODUCTION 

Quantum transport phenomena that involve the car- 
rier's spin degree of freedom have attracted a lot of at- 
tention during the past decade. 0,0 A number of stud- 
ies have analyzed the spin polarized transport in ballis- 
tic regime and reported intriguing phenomena, 0,3 and 
have opened up the possibility of new spintronic devices. 
It is also interesting to study the spin-dependent trans- 
port in diffusive and chaotic regimes because the inter- 
ference of coherent electron waves shows a number of 
characteristic effects when the spin degree of freedom is 
taken into consideration. One of the special character- 
istics of diffusive and chaotic systems is the fluctuating 
nature of transport coefficients such as conductance. Q It 
is well understood that such fluctuations do not depend 
on the details of the sample parameters, but depend only 
on the symmetry of the system. 0, |^ The two relevant 
symmetries here are time-reversal symmetry (TRS) and 
spin-rotation symmetry (SRS). TRS is broken by applied 
magnetic fields or by magnetic scattering due to magnetic 
impurities or magnetic domain walls. If TRS is broken, 
systems are classified as unitary, regardless of whether 
or not SRS is broken. The spin-orbit interaction breaks 
SRS but preserves TRS, and in this case the systems are 
classified as symplectic. 

When the spin degree of freedom is taken into ac- 
count, the description of the conductance fluctuations 
becomes more complex. Conductance fluctuations of 2- 
dimensional systems which are coupled to the Ising spin 
glass is reported. 7] It has been reported that the reduc- 
tion of variance of conductance takes place due to the 
Zeeman splitting in the sample region. H 1 [13 Ait- 
shuler and Shklovskii had shown that the variance of the 
conductance is described by@ 



Recent works have pointed out that the universal 
conductance fluctuations (UCF) in a chaotic quantum 
dot in t he p resence of spin-orbit scattering shows new 
features. [Hi [HE! They have shown that the UCF in 
the presence of spin dependent scattering is interesting 
not only from the theoretical point of view but also from 
the experimental view point. 

The transport properties of mesoscopic systems de- 
pend not only on the sample but also on the states in 
the leads through which currents flow into and out of 
the sample and through which voltages are measured. |l4j 
How the transport properties and the universality class 
are changed by modulating lead states is a very interest- 
ing question, especially when the spin degree of freedom 
plays a role. 

In this paper, we investigate the influence of the spin 
scattering on transport properties in disordered systems. 
We consider magnetic impurities in sample region and the 
Zeeman field in the lead. Three types of magnetic impu- 
rities are considered. We call these Ising, XY and Heisen- 
berg, depending on the number of spin components. In 
order to calculate two terminal conductance, we employ 
the transfer matrix methodjl^l that is extended to in- 
clude the spin degree of freedom. Magnetic impurities 
remove the spin degeneracy and break TRS in certain 
cases. 

We find that the variance of the conductance is halved 
as the number of spin components of the magnetic im- 
purities increases. When the Zeeman field is applied in 
a lead, a further reduction of the variance is observed. 
In order to observe the crossover of the universality class 
with the increase of the Zeeman field in the lead, we study 
the level spacing distribution of transmission eigenvalues. 
Part of this work has been presented in the international 
conference, "Localisation 2002". OE3 



(i) 



where bd is a dimension-dependent factor that is of the 
order of unity. (3 is equal to 1,2 and 4 for orthogonal, 
unitary and symplectic systems, respectively. The quan- 
tity K is equal to the number of noninteracting series of 
energy levels with s-fold spin degeneracy. 



II. MAGNETIC SCATTERING 

We consider a two dimensional (2D) system connected 
to two electrodes. The 2D system is constructed in the x- 
and y- directions and the current flows in the x direction. 
There is an exchange interaction between the electron 
spin and the static local spins in the system. The one- 
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FIG. 1: Variance of conductance in the presence of magnetic 
impurities. System size is 30 x 30. W is set to be 4.0 and 
J = 0.4. 



FIG. 2: Variance of conductance when the sample contains 
both the magnetic field and magnetic impurities. Parameters 
are the same as in Fig. 1. 



electron Hamiltonian is 

H = H + H', (2) 

H Q = ^2Wicf a C ita - V ic,jcr'4,<T C j,<T'i ( 3 ) 

h' = -j J2 c L°v-' c ^' • ( 4 ) 

where c ia .(ci >cr ) denotes the creation (annihilation) oper- 
ator of electron at the site i with spin <r(= ±1) on the 
2D square lattice. On-site energy W, denotes the non- 
magnetic random potential distributed independently 
and uniformly in the range [—W/2,W/2]. The hopping 
is restricted to nearest neighbors. 

We investigate how the variance of the conductance 
distribution of Ho changes due to the presence of H' , 
The variance of conductance distribution for Hq is de- 
termined by the symmetry, which is controlled by the 
hopping term, Via-ja' ■ If it is set to 2 x 2 unit matrix, 
H belongs to the orthogonal class. H belongs to the 
unitary class when 

_ / exp(i</>) \ , , 

V ^+y -\ exp(i0) J ' [b) 

where <fi is a Peierls phase which is distributed uniformly 
in the range [0, 2 tt). To realize systems belonging to the 
symplectic class, [l8| we- set the hopping 

/ cos6> sin6> \ , . 

V ^={-sm8 cose ) ^ 

and 

_ / cos 9 i sin 6 \ , , 

Vi ' t+ y ~ yisme cos6> ) ■ [ '> 



where the parameter 8 denotes the strength of the spin- 
orbit interaction. 

The additional term H' is the spin scattering term, a 
is the Pauli spin matrix and Si is the static local spin. 
We consider three types of magnetic impurities: Ising, 
XY and Heisenberg where Si = (0, 0, S z ), (S x , S y ,0) and 
(S x , S y ,S z ), respectively. The absolute value of Si is set 
to unity, and their components are distributed randomly 
according to the uniform distribution on n-dimensional 
sphere, n(=l,2 or 3) being the number of spin compo- 
nents. 

In order to calculate the conductance, we employ the 
transfer matrix method 15] and extend it to include the 
spin degree of freedom. The conductance G is given by 
the Landauer formula |l9ll2fj| | as 

G = ^tr(ttt) = f!^ Tl ( 8 ) 

i 

where t is the transmission matrix including the spin de- 
gree of freedom and Tj is the transmission eigenvalue. In 
the present simulation, the system size is 30 x 30 in units 
of the lattice spacing 

In Fig.l we show the variance of conductance in the 
presence of magnetic impurities where Hq belongs to the 
orthogonal class. W is set to be 4.0 and J = 0.4. 10 4 
ensemble averages are taken for each data. When J = 0, 
there is no magnetic scattering and the variance is close 
to that expected for the universal conductance fluctu- 
ation of orthogonal systems. Q For convenience, let us 
denote the variance of Tj for the conventional orthog- 
onal class as V or th and G/ (e 2 /h) — G, so that in the case 
of J = 0, VarG = V OIt h- Votth. is 1/2 for the chaotic cav- 
ity, 8/15 for the quasi-lD wire, 0.74 for 2D, and 1.2 for 
3D systems, respectively. 1(| 
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FIG. 3: Variance of conductance when the sample contains 
both the spin-orbit interaction and magnetic impurities. Pa- 
rameters are the same as in Fig. 1. 



FIG. 4: The effect of the Zeeman coupling in a lead on the 
variance. Ho belongs to the orthogonal class, W = 3.0 and 
Z = 3.0 



In the presence of the Ising type impurities, no spin 
flip processes occur as in the conventional non-magnetic 
system. However, the up spin and down spin electrons 
are described by different wave functions, because the 
exchange field of magnetic impurities lifts the spin de- 
generacy. With sufficiently large J, these wave functions 
are no longer correlated, and the variance is represented 
by 

VarG = VarG T + VarG x , (9) 

where Gjm is the conductance through the up (down) 
spin channel. Since both of the variances Gj and G x are 
Vorth/4) the sum G becomes V ort h/2 

While Ising type impurities do not rotate the spin di- 
rection, the spin flip process occurs in XY type impu- 
rities. Then the variance is simply given by VarG = 
Vorth/4, since the factor of 4 coming from the spin degen- 
eracy in the conventional non-magnetic orthogonal class 
is missing. Though the Hamiltonian is complex due to 
a y , the statistics of the transmission eigenvalues as well as 
the energy level statistics are that of the orthogonal class. 
This can be seen if we define the time-reversal operator 
T as T = i<J x K, where K denotes the complex conjuga- 
tion. This operator is anti-unitary, satisfies T 2 = 1 and 
commutes with the Hamiltonian of the system including 
XY type impurities. 

For Heisenberg and XY type impurities, spin flips oc- 
cur. However, when H' includes the Heisenberg type 
scatterers, the Hamiltonian no longer commutes with the 
time-reversal operator, and the system is classified into 
the unitary class. Therefore, the variance is further re- 
duced by a factor of 2. 

From these results, the variance of conductance in the 



presence of impurities is given by 

VarG=^, (10) 

We have numerically inverstigated a square 2D system, 
but the argument above is general, and this relation 
should be valid in other dimensions. [16| 

Figure 2 shows the change of the variance when Hq 
is in the unitary class. Parameters are the same as in 
Fig. 1. We consider the random magnetic field and the 
phase in the transfer integral is given by Eq. (2). Without 
magnetic impurities, the system belongs to the unitary 
class with spin degeneracy and the variance is Vorth/2. 
Ising type magnetic impurities removes the spin degen- 
eracy and the variance becomes Vorth/4. Spin flips occur 
due to the XY type magnetic impurities and the vari- 
ance becomes half of the Ising case. As shown in Fig. 2, 
the variance of the conductance when H' includes the 
Heisenberg type impurities is the same as in the case of 
the XY type, because the Hamiltonian of neither system 
commutes with the time reversal operator. 

The reduction of the variance is also obtained when Hq 
includes the spin-orbit interaction and the Hamiltonian 
belongs to the symplectic class (Fig. 3). Once the mag- 
netic impurities are included, irrespectively of the type 
of the magnetic scatterers, the variance becomes Vorth/8 
from F or th/4. 



III. EFFECT OF LEADS 

We then address the question of whether or not trans- 
port properties and the universality class can change as a 
result of an asymmetric spin population in the leads, i.e., 
when the number of up and down spin channels in the 
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FIG. 5: Spacing distribution of r for the sample with the 
spin-orbit interaction. We set W = 2.0 and Z — 3.0. The 
distribution without Zeeman splitting in a lead (o) fits to the 
form of GSE while that with Zeeman (•) fits to GUE. 




FIG. 6: Schematic diagram of the variance of the conduc- 
tance. The horizontal arrow indicates the increase of the spin 
component of the magnetic scatterers, while the slanted ar- 
row indicates the addition of the Zeeman field. The Zeeman 
field is assumed to contain a component different from that 
of the magnetic scatterers. 



leads becomes asymmetric as a result of Zeeman split- 
ting. We show that even if the sample is unchanged, the 
universality class can change under certain conditions. 

We consider the Zeeman field in one of the leads, set- 
ting the other lead Zeeman field free. The transverse 
energy of the channel (i, a) in the lead ef is given by 

ITT 

e f = -2cos(-^-)-Z<7 (i=l,2,-.. ,L) (11) 

where L denotes the number of sites in the transverse 
direction (y-direction), and Z denotes the strength of 
Zeeman splitting and er(= ±1) is the spin index. The 
channel (i, a) is a propagating mode if \Ep — ef\ < 2.0, 
Ep being the Fermi energy. For example, when we set 
Ep = —1.1 and Z = 1.0, the number of up (down) spin 
channels becomes 27 (15) for the sample of width 30 sites. 

Figure 4 shows the variance in the presence of Ising 
and XY type impurities in a system including Zeeman 
splitting in a lead. In this simulation, we set W = 3.0, 
9 = 7r/4 and Ep = —1.1. The system size is again set to 
be 30 x 30 in units of the lattice spacing. The population 
of up and down spins in one lead is always set to be 
symmetric [Z = 0), and that in the other lead is varied 
(Z = 0,3.0). 

From this figure, we observe that the variance becomes 
almost half due to the Zeeman splitting in the leads, even 
if the sample region is not changed at all. The effect 
of Zeeman splitting is similar to that of increasing the 
number of spin component of the magnetic impurities. It 
should be noted that if the direction of the Zeeman field 
is the same as one of the spin components of the scatters, 
the reduction in not observed. For example, in the case of 
Ising type scatterers, applying the Zeeman field in the z- 
direction does not change the variance while the halving 
of the variance is observed if we apply the Zeeman field 



in x- or y-direction. The same is true for the case of the 
XY type scatters. We need to apply the Zeeman field in 
the ^-direction to observe the halving of the variance. 

We then show that the Zeeman field in the lead changes 
the universality class of the system. To detect the 
crossover of the universality classes, we investigate the 
spacing distribution P(s) where s is the interval between 
neighboring transmission eigenvalues r's. An ensemble 
of about 10 6 samples is simulated to get good statistics. 

Figure 5 shows the spacing distribution P(s) of the 
transmission eigenvalue r for the sample with spin-orbit 
interaction. In the absence of Zeeman splitting in the 
leads (o), P(s) is close to the Wigner surmise for the 
Gaussian symplectic ensemble (GSE). On the other hand, 
with Zeeman splitting in the leads (•), P(s) is close to 
that for the Gaussian unitary ensemble(GUE). The rea- 
son for this crossover is that the asymmetry of up and 
down spins in the lead destroys the self-dual property of 
the scattering matrix. To be specific, we consider the 
(i, j) component of the S-matrix Sij which describes the 
transmission from (i |,i J,) to (j \,j j). The self-duality 
requires 

M^.M-c?) <12) 

which means that the transmission probability from i "f 
to j | is the same as that from j | to i J,. This symmetry 
is broken when we apply the Zeeman field in the lead, 
which lifts the degeneracy of i 1 and i |. 
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TABLE I: Variance of the system in the presence of magnetic impurities and the Zeeman field. 



Universality class of Ho 


Additional spin scattering 


Zeeman field in a lead 


Variance/Forth 


Orthogonal 





no 


1 







yes 


1/2 




Ising 


no 


1/2 




Ising 


yes 


1/4 




XY 


no 


1/4 




XY 


yes 


1/8 




Heisenberg 


irrelevant 


1/8 


Unitary 





no 


1/2 







yes 


1/4 




Ising 


no 


1/4 




Ising 


yes 


1/8 




XY 


irrelevant 


1/8 




Heisenberg 


irrelevant 


1/8 


Symplectic 





no 


1/4 







yes 


1/8 




Ising 


irrelevant 


1/8 




XY 


irrelevant 


1/8 




Heisenberg 


irrelevant 


1/8 



IV. SUMMARY 

We have studied the effect of spin scattering in dis- 
ordered systems on the fluctuating nature of the con- 
ductance. We have considered magnetic impurities in a 
sample and calculated transport properties. Our results 
show that the variance of conductance is halved as the 
number of spin components of the magnetic impurities 
increases. Halving of the variance of the conductance is 
also obtained when the sample includes magnetic field or 
spin-orbit interaction. 

We have also investigated the effect of the Zeeman 
splitting in a lead. Halving of the variance of conduc- 
tance is obtained when the direction of the Zeeman field 
contains the component different from the component (s) 
of the magnetic scatterers. This behavior is reminiscent 
of the change of the variance in the superconducting- 
normal junction. 0, l2ll |2^| Analyzing the transmission 
eigenvalues, the universality class has been shown to be 
changed by the Zeeman field in the lead. The results are 
summarized in Table I and schematically shown in Fig. 6. 

Before concluding, we relate our results with that 
of Aleiner and Fal'ko.[l3j They have obtained for the 
chaotic system 



where s — 1,2 indicates the Kramers degeneracy and 
E = 2 if the spin flip process is present and E = 1 other- 
wise. (3 =1,2 or 4 is determined by the universality class. 
Setting V^ r th = 1/2 we recover their interesting result. 
For example, when Hq is classified into the orthogonal 
class, and H' includes the XY type magnetic impurities, 
/3 = 1, s = 1 and E = 2, which gives VarG = | = Krth/4. 
Therefore, the present results are the extension of Aleiner 
and Fal'ko [l3| to higher dimensions and to the inclusion 
of the effect of the Zeeman splitting in the lead. 



VarG 



4/3E' 



(13) 
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